Infrared dynamic polarizability of HD + rovibrational states 



J.C.J. Koelemeij* f 



O 

(N 

C 

Oh 
i 

a 

o 
-i— > 

O 



A calculation of dynamic polarizabilities of rovibrational states with vibrational quantum number v = — 7 and rotational 
quantum number J = 0, 1 in the lscr g ground-state potential of HD + is presented. Polarizability contributions by transitions 
involving other lscr ff rovibrational states are explicitly calculated, whereas contributions by electronic transitions are treated 
quasi-statically and partially derived from existing data [R.E. Moss and L. Valenzano, Molec. Phys., 2002, 100, 1527]. Our 
model is valid for wavelengths > 4 fim and is used to to assess level shifts due to the blackbody radiation (BBR) electric field 
encountered in experimental high-resolution laser spectroscopy of trapped HD + ions. Polarizabilities of 1sct s rovibrational states 
obtained here agree with available existing accurate ab initio results. It is shown that the Stark effect due to BBR is dynamic 
and cannot be treated quasi-statically, as is often done in the case of atomic ions. Furthermore it is pointed out that the dynamic 
Stark shifts have tensorial character and depend strongly on the polarization state of the electric field. Numerical results of BBR- 
induced Stark shifts are presented, showing that Lamb-Dicke spectroscopy of narrow vibrational optical lines (~ 10 Hz natural 
linewidth) in HD + will become affected by BBR shifts only at the 10 -16 level. 

1 Introduction 
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The molecular hydrogen ion (H^) and its isotopomers 
(HD + , D^, HT + , etc.) are the simplest naturally occurring 
molecules. As such they are amenable to high-accuracy ab ini- 
tio level structure calculations, which are currently approach- 
ing 0. 1 ppb for rovibrational levels in the electronic ground 
potential 1 . The inclusion of high-order QED terms in these 
calculations makes molecular hydrogen ions an attractive sub- 
ject for experiments aimed at comparison with theory and tests 
of QED. With rovibrational states having lifetimes exceeding 
10 ms it has long been recognized that optical (infrared) spec- 
troscopy could provide accurate experimental input, and sev- 
eral experimental studies were undertaken 2 3 or are currently 
in progress 4 . The highest accuracy that has hitherto been 
achieved is 2 ppb for a Doppler-broadened vibrational over- 
tone transition at 1 .4 /im in trapped HD + molecular ions, sym- 
pathetically cooled to 50 mK 3 . By comparison, the highest ac- 
curacy achieved in laser spectroscopy of laser-cooled atomic 
ions, tightly confined in the optical Lamb-Dicke regime, is 
~ 1 x 10" 17 in the case of the A1+ optical clock at NIST 
Boulder, USA. 5 . The Al + optical clock employs quantum- 
logic spectroscopy (QLS) which utilizes entangled quantum 
states of two trapped ions, one of which is used for (ground- 
state) laser cooling and efficient state detection, whereas the 
other ion contains the transitions of spectroscopic interest 6 . It 
has been pointed that Doppler-free spectroscopy may be per- 
formed on HD + as well 3 , and also that QLS may be used for 
spectroscopy of molecular ions 7 . 
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Fig. 1 (Color online)(a) Potential energy curves of the ls<r 9 and 
2pa u electronic sates. Indicated energy values are binding energies 
of the molecule. Shown also are radial nuclear (vibrational) 
wavefunctions, Xv(R), for v — 0, 4 as well as one dissociating 
nuclear wavefunction in the 2pa u state. The red arrow represents a 
purely rovibrational transition within lscr fl ; the blue arrow 
exemplifies a transition between different electronic states, (b) 
Dipole moment function Di (R) used in the calculation of radial 
dipole matrix elements (solid curve), shown together with the 
approximate function used by Colbourn and Bunker 8 (dashed line) 
and the fully g/u symmetry-broken dipole moment function, valid at 
long internuclear range (dot-dashed line). 



Accurate results of laser spectroscopy of HD + are of inter- 
est for the determination of the value of the proton-electron 
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mass ratio, m p /m e 2 , and for the search for a variation of 
m p /m e with time 9 . The former may be achieved by combin- 
ing ab initio theoretical results with results from spectroscopy 
at an accuracy level of ~ 10~ 10 ; for the latter spectroscopic 
results with an accuracy of ~ 10~ 15 are required to improve 
on the current most stringent bounds 10,11 . In both cases spec- 
troscopy of optical transitions is faced with level shifts due to 
magnetic and electric fields and, to a lesser extent, shifts due to 
collisions and relativistic effects. The Zeeman effect of HD + 
was recently considered by Bakalov et ai, and level shifts to 
second order in the magnetic field were given for a large set of 
rovibrational states 12 ' 13 . Static polarizabilities of vibrational 
states with rotational quantum number J = 0, 1 were calcu- 
lated and reported by several authors 14 ' 15 , while dynamic po- 
larizabilities of HD + vibrational states with J = were eval- 
uated for a discrete set of two-photon transition wavelengths 
in the 1-18 ^m wavelength range 15 . However, to our knowl- 
edge, no results on dynamic polarizabilities of HD + for vi- 
brational states with J > are available in literature. Polar- 
izabilities of such states for a wide range of infrared wave- 
lengths are required for the calculation of differential Stark 
shifts due to blackbody radiation (BBR). Moreover, since the 
BBR spectrum encompasses several rovibrational transitions 
of the HD + ion, it is expected that the quasi-static treatment of 
BBR-induced Stark shifts as often done in the case of atomic 
ion species is not valid for HD + . Rather, the case of HD + 
will be analogous to that of neutral polar molecules, for which 
BBR-induced Stark shifts were evaluated using dynamic po- 
larizabilities 16 . 

This Article addresses the (BBR-induced) dynamic Stark 
effect of HD + and is organized as follows. In Sec. 2 we 
present our model to calculate dynamic polarizabilities and 
BBR-induced Stark shifts, followed by a discussion of the re- 
sults for several rovibrational states in the 1sct s ground-state 
potential of HD + in Sec. 3. Conclusions are presented in 
Sec. 4. Throughout this Article, the terms 'Stark effect' and 
'polarizability' will be used interchangeably, and SI units will 
be used. 

2 Theory 

Figure 1(a) shows a partial energy level diagram of the HD+ 
molecular ion including the electronic ground-state poten- 
tial, lsa g , and the first electronically excited potential, 2per„. 
Note that in HD+ the g/u symmetry quantum labels are only 
approximately good quantum labels as the nonidentical nu- 
clei introduce g/u symmetry breaking at large internuclear 
range. The potential energy curves shown are interpolations 
of data published by Esry and Sadeghpour 17 , who present the 
potential energy as the sum of a nonrelativistic, fully adia- 
batic curve, and a diagonal nonadiabatic correction. We use 
these curves to obtain (real-valued) radial wavefunctions of 



nuclear motion, Xvj(R), by numerical solution of the radial 
Schrodinger equation including the centrifugal term due to the 
molecular rotation: 
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Xvj(R) = e v jXvj{R), 

(1) 



where R denotes the internuclear separation, fi stands for the 
nuclear reduced mass of the molecule, v labels the vibrational 
state, J is the rotational angular momentum of the molecule, 
and E v j is the rovibrational energy. Vi(R) are the potential 
energy curves for the i — lscr g , 2pa u states taken from Esry 
and Sadeghpour 17 , who also provide dipole moment functions 
Di(R) and Di 2 (R). These correspond to the dipole moment 
of the 1 sa g state and the dipole moment of electronic transi- 
tions between lser g and 2pa u , respectively. 

The dynamic polarizability corresponds to the ability of the 
HD + molecule to deform under the influence of an oscillating 
electric field, and depends on the strengths and frequencies 
of many electric dipole transitions in both the nuclear and the 
electronic degrees of freedom. Laser spectroscopy on HD + 
is typically performed on transitions between low-lying rovi- 
brational levels in the lscr g state, and it is the dynamic polar- 
izability of these levels that we will focus on here. The dy- 
namic polarizability, a(oj), is defined as follows. A quantum 
state with quantum numbers (v, J, M) (with M corresponding 
to the projection of J on the space-fixed z-axis) and energy 
E v .jm will undergo an energy shift AE due to the interaction 
with a monochromatic electric field with amplitude £, polar- 
ization state q, and angular frequency co equal to 



AE = - 



1 



l vJM 



(u)£ 2 



(2) 



In the remainder of this Article, the polarization state q 6 
(—1,0,1) will be taken q = (i.e. linear polarization parallel 
to the space-fixed z-axis) and we will omit the label q alto- 
gether; see Sec. 2. 1 . The dynamic polarizability can be written 
as a v j M (w ) = < j M (w ) + a% JM (w ) , where a r v v JM {u) stands 
for the contribution by transitions coupling to other lser g rovi- 
brational states, and a.% JM (u>) accounts for the contributions 
by all transitions connecting to electronically excited states. 
To simplify the calculation, we will restrict ourselves to the 
two strongest sets of transitions from the rovibrational states of 
interest. These are (1) purely rovibrational transitions within 
the electronic ground state lser g , and (2) electronic dipole tran- 
sitions to dissociating states in 2pa u . 

2.1 Dynamic polarizability due to rovibrational transi- 
tions in lsag 

As a starting point we will use the energy shift derived using 
time-dependent second-order perturbation theory, [Eq. (7.73)] 
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Table 1 Static polarizabilities (in units of 47reoOo) for vibrational 
states with J = 0. Total polarizabilities a vJM (0) were taken from 
Moss and Valenzano 14 . Individual rovibrational and electronic 
contributions alj M (0) and o^ JM (0), respectively, are also 
specified. Entries in the rightmost column are obtained from those 



in the other columns as a v jM 


(0) - aZ M (0) 




v cwm(0) (Ref. [14]) 


o«jm(0) (this work) 


"Sjm(0) 


395.306 


392.2 
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1 462.65 
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3.8 


2 540.69 
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4.5 


3 631.4 
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5.5 


4 737.3 


730.7 


6.6 
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6 1008 


998.6 
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7 1184 
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Table 2 Comparison of purely electronic static polarizabilities 
o^ JM (0)(in units of 47reoao) f° r vibrational states with J = of 
HD + with accurate values for vibrational states with J = of 
and Dj, calculated by Hilico et al. 18 Note that for and the 
static polarizability stems from electronic transitions only 
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16.290 999 14 
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9.162209 589 



in the textbook by Sobelman , 

(3) 

where v,J,M and v',J', M' are the quantum numbers of the 
initial and final states, respectively. Here we assume that the 
states (v, J, M) are degenerate in the quantum number M. It 
is furthermore important to note the role of the sign in the 
definition of to vv ijjr. 

Ww'JJ' — {E v 'ji - E vJ )/h. (4) 

Hence, uj vv >.j.j' > for transitions to more highly-excited 
states and uj vv >ll' < for transitions to lower states. For 
purely rovibrational transitions, the squared dipole transition 



matrix element \D vv >jj'MM'q\ 2 reduces to 20 

\D vv ..j.r M M- q \ 2 = \{JM\V*_ q0 (cj E )\J'M')\^ 2 vvlJJ , 

= (2J + D(2J' + 1)(^ J ') 2 

( J 1 J' V 2 

X \ -M -q M'J Vvv'JJ', (5) 

with 

poo 

Hlv>jj> = \ xv.r{R)D 1 {R)xv,j{R)dR\ 2 . (6) 
Jo 

The dipole matrix element li vv <jj> is a vector oriented along 
the internuclear axis of the HD + molecule. Therefore, in or- 
der to evaluate the matrix elements, li vv <jj> needs to be trans- 
formed from the molecule-fixed to the space-fixed frame by 
rotation about the set of Euler angles oje, which is imple- 
mented through the rotation operator X>* g0 (wE) in the first 
factor in Eq. (5). In arriving at the second line of Eq. (5) we 
use the fact that for states with A = (like for ls(T g , while 
ignoring the spins of the proton, deuteron and electron) the 
projection of J on the internuclear axis is zero. As stated in 
Sec. 2, we will consider the case q — only. 

The squared matrix elements ii 2 v ,jj, are readily evalu- 
ated using the numerical expressions for wavefunctions and 
dipole moment functions introduced above. The expression 
for 

a JjjM( w ) * s obtained after inserting Eqs. (4) and (5) into 
Eq. (3), followed by equating Eq. (3) to Eq. (2) and solving for 
a vjM( u} ) (momentarily assuming that a^ JM (uj) = 0). As we 
here focus on low-lying vibrational levels and dipole transi- 
tions only, we will truncate the summation in Eq. (3) to v = 9, 
and also ignore the contribution by purely rovibrational tran- 
sitions to continuum states above the lscr ff dissociation limit. 
This is justified as the line strength of vibrational overtones 
decreases rapidly with increasing order of the overtone. The 
summation is furthermore limited to terms obeying the selec- 
tion rule J' = J ± 1. 

2.2 Polarizability due to electronic transitions 

For static electric fields (oj — > 0), it is known that a r v v ]M {u) 3> 
a v.jM ( w ) 14 ■ This may not necessarily be the case for infrared 
frequencies, for which a™ JM (uj) is expected to be smaller as 
spectrally nearby vibrational overtones are generally weak, 
whereas the detuning from strong rotational transitions and 
fundamental vibrations is large. Thus, there may be spectral 
regions where a e vJM {uj) becomes comparable in magnitude to 
a£j M (w). However, transitions from low-lying lser g rovibra- 
tional states to 2pcr„ states are located in the ultraviolet (UV) 
or even in the vacuum-ultraviolet (VUV) spectral range. Since 
the frequencies present in the T = 300 K BBR spectrum are 
in the infrared (peak emission wavelength ~ 10 ^m), it seems 
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Table 3 Static polarizabilities (in units of 47reoao) for vibrational states with J = 1. Total polarizabilities owm(O) were taken from Moss 
and Valenzano 14 . Individual rovibrational and electronic contributions 

d*v v jm (0) «md ck!^j^(0), respectively, are also specified. For each DA 
value, entries in the rightmost column are obtained from those in the other columns as cwm(0) — a£j M (0) 
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-274.0 


5.1 


141.50 
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220.3 
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-497.57 
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263.22 
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5.9 
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-580.99 


-596.7 


16 
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301.0 


7.1 
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-679.82 


-699.9 


20 


361.66 


353.0 


8.7 



justified to regard the BBR electric field as static where it con- 
cerns al JM (uj). In Sec. 3.3.1 it will be further justified that for 
this reason, a£ JM (0) is a good approximation to a% JM {uj). 

Rather than deriving the static polarizability a^ JM (0) 
from second-order perturbation theory, we extract its values 
from previously published and accurate static polarizabilities 
a vjM (0)' obtained by a full nonadiabatic calculation by Moss 
and Valenzano 14 , as follows. From each of the total static po- 
larizabilities C( v jm(0) tabulated by Moss and Valenzano, we 
subtract our value for a r v V ]M (0) calculated using the procedure 
described in Sec. 2.1 to obtain a e uJM (Q). 

3 Results and discussion 

3.1 Rovibrational wavefunctions and dipole matrix ele- 
ments 

Before discussing the results of our method to obtain 
oi v .jM{u), it will be worthwhile to investigate the accuracy 
of the wavefunctions Xvj(R) and energy levels E v j obtained 
from Eq. (1), as well as the accuracy of the radial dipole matrix 
elements li vv <jj> calculated using Eq. (6). From comparisons 
with more accurate nonrelativistic level calculations for HD + 9 
the inaccuracy of the energies E v j calculated here is found to 
be a few parts in 10 5 (or less than 0.5 cm^ 1 ), in correspon- 
dence with the accuracy specified by Esry and Sadeghpour 17 . 
The accuracy of the energy levels also gives an indication of 
the accuracy of the wavefunctions Xvj(R)- 

In order to check the accuracy of the radial matrix elements 
Hw'jj', a comparison can be made with values calculated by 
Colbourn and Bunker 8 . Here it is important to note, how- 
ever, that Colbourn and Bunker ignore effects of g/u symmetry 
breaking by using a dipole moment function Dqb (R) ~ ei?/6 
(with e the electron charge) This functional form is valid 



§This expression follows from evaluating the HD + ls<r 9 dipole moment with 
respect to the center of mass at the equilibrium internuclear separation, for 
which the electron on average sits halfway the two nuclei. 



at short internuclear range, where effects of g/u symmetry 
breaking are small. However, for large internuclear separa- 
tion in the lser g state of HD + , the electron sits primarily at the 
deuteron, which leads to a dipole moment function varying 
for large R as w (2/3)e_R. The function Di{R) provided by 
Esry and Sadeghpour includes effects of g/u symmetry break- 
ing, as illustrated in Fig. 1(b). To compare with the results by 
Colbourn and Bunker, we first use our Xvj{R) with Dqb(R) 
to obtain matrix elements P^ijji. We find agreement at the 
level of a few times 1CT 5 , consistent with the accuracy of both 
our wavefunctions Xvj(R) and those used by Colbourn and 
Bunker, which produce energy levels with similar accuracy. A 
second calculation using D\{R) instead of Dcb(R) leads to 
radial matrix elements differing from those by Colbourn and 
Bunker at the level of 2 x 10~ 3 for transitions v' = 1 — v = 0, 
and 4 x 10~ 3 for v' = 5 — v = 4. This difference we attribute 
to the inclusion of g/u symmetry-breaking effects in D\{R), 
and may be considered an improvement over the values by 
Colbourn and Bunker. We put an conservative error margin 
of 25% on this difference, thereby placing an upper bound of 
1 x 10~ 3 on the accuracy of the matrix elements \i vv 'jj> . 

3.2 Static polarizability results 
3.2.1 Accuracy of 

a »jif(^)' The results of Sec. 2.1 en- 
able us to calculate dynamic polarizabilities a™ JM {Lo). To as- 
sess the accuracy of these calculations, we have checked the 
dependence of the static polarizability a™ JM (0) on the accu- 
racy of both the energy levels and the radial matrix elements 
used. Computing a r v v JM (Q) once with the eigenvalues E v j 
of Eq. (1), and once with accurate energy levels published 
by Moss 21 (accuracy better than 0.001 cm -1 ), we find that 
a£j M (0) varies by ~ 1 x 10~ 4 . A similar check is done by 
using values \/i vv >jj' | 2 computed using Dcb(R) and Di(R), 
respectively. The effect of the improved values on a r v v JM (0) 
is a few times 10~ 3 . Placing again a conservative bound of 
25% on the accuracy of this improvement, the accuracy of our 



4 



1-9 



03 

o 

to 



500 r 



50 ; 



10" 

& 5 





(7,0) -(9,1) (0,0) -(1,1) 




A 


■• — *- 


" • i 

1 
1 








. . . 




(1,0) -(0,1) ; (7,0) -(6,1) 
(1.0) -(2,1) ! (7,0) -(8,1) 



6 7 

wavelength / |^m 



10 



11 



Fig. 2 (Color online) Dynamic polarizabilities versus wavelength (A = 2hc/oj) for states with J = and, from bottom to top, v = 0, 1, 7, 
respectively. The curves were produced using Eq. (7). For each curve, dashed segments correspond to negative values, while solid segments 
correspond to positive values. Each 'dip' or 'peak' in a curve corresponds to a zero crossing of the polarizability. Vertical dashed lines indicate 
the position of rovibrational transitions (v,J) — (V, J') coupling to J' = 1 states. The dots show the more accurate values calculated for 
specific wavelengths by Karr et at 15 , which agree with our result to within 3%. 



value of a£j M (0) is found to be < 1 x 10~ 3 . We also moni- 
tored the effect of the truncation of Eq. (3) to v' = 9. This has 
no noticeable effect at the 1 x 10 -3 level for states v < 7. 



owm(w) = a r v v JM (uj) + o% JM (u) by the expression 



(7) 



3.2.2 Accuracy of ot% JM {G). As described in Sec. 2.2, 
the values ct^ v JM (0) may be combined with previously pub- 
lished values qwm(O) to extract o^jm(O)- Thus-found val- 
ues of ct e vJM (0) are presented in Tables 1 and 3. We find that 
a vjM^) contributes to <Xu,/m(0) at the 1% level. Given the 
< 1 x 10 -3 accuracy of our results for a r v v j M (0), we are lead 
to believe that the values of aJj M (0) inferred here are accu- 
rate to within 10%. 

It is furthermore interesting to compare the values of 
aJj M (0) obtained here with static polarizabilities of the iso- 
topomers and D^, which were calculated with high accu- 
racy for vibrational states with J = by Hilico et al. 1S . In 
Table 2 it can be seen that for each vibrational state, the HD + 
value lies in between the values for and D J . This is ex- 
plained by the fact that the energy of a given vibrational state 
scales as J\j 'fi, with /i the reduced nuclear mass of the iso- 
topomer. Thus, for large reduced mass, vibrational levels are 
more deeply bound and therefore exhibit a smaller static polar- 
izability. As the variation of binding energy is small compared 
to the typical energies of transitions to 2p<7 u states, the mass 
scaling of the polarizability is approximately linear, and the 
value for HD + should be located halfway the values for 
and T)^ as m Table 2. 

3.3 Dynamic polarizability results 

3.3.1 Accuracy of the approximation. As discussed 
in Sec. 2.2, we will approximate the dynamic polarizability 



For the infrared spectral range of interest here (A > 4 /jm) 
we believe that by approximating o^j M (o;) by a e vJM (Q) we 
systematically underestimate the magnitude of the shift due 
to al JM (uj) alone by less than 10% (details of this estimate 
are postponed to the Appendix). This is comparable to the 
uncertainty of the values al JAI (0) reported in Tables 1 and 
3. In order to verify the accuracy, we compare the result of 
Eq. (7) with the more accurate values calculated by Karr et al. 
for a discrete set of wavelengths for states with J = (Fig. 2). 
The results of the two methods are found to agree within 1 % 
for v — and within 3% for v = 7. As the comparison is made 
for relatively short wavelengths, for which the polarizability 
stems almost entirely from al JM (uj), the level of agreement 
is consistent with the estimated error of < 10% in the value of 
<jm(0)- 

The result for a^j M (0) obtained here is more useful than 
one would expect on the basis of its error margin for two rea- 
sons. First, for dynamic Stark shifts due to BBR (found by 
integrating the dynamic Stark shift over the BBR electric field 
spectral density; see Eq. (9) and the Appendix), we estimate 
the error introduced by the quasi-static approximation to be 
even smaller, < 3%. Second, for spectroscopy one is primar- 
ily concerned with differential level shifts, for which the sys- 
tematic errors in q: v jm(uj) will partially cancel. 

3.3.2 Dependence on \M | and polarization state. It was 

mentioned in Sec. 2.1 that Eq. (7) tacitly assumes linearly po- 
larized electric fields. For obtaining the shift due to unpolar- 
ized, incoherent BBR, it is necessary to average over the three 
independent polarization states q = —1,0,1. It may be shown 
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Fig. 3 (Color online) (a) Dynamic polarizabilities versus wavelength (A = 2-kc/lj) for various states with v — 0, J — 1, computed using 
Eq. (7). For each curve, dashed segments correspond to negative values, while solid segments correspond to positive values. Each 'dip' or 
'peak' in a curve corresponds to a zero crossing of the polarizability. Vertical dashed lines indicate the position of rovibrational transitions 
(v, J) — (v', J') coupling to J' — 0, 2 states. The curves show marked tensorial differences between different A/-states for polarized electric 
fields. It is also seen that for shorter wavelengths the contribution by rovibrational transitions becomes less significant, and that the electronic 
contribution becomes dominant instead. Furthermore, the magnitude of the average polarizability is seen to decrease towards longer 
wavelengths, which can be inteipreted as an geometric averaging effect of the molecular rotation, (b) Mean-square electric field spectral 
density of the BBR at T — 300 K. The BBR spectrum encompasses several rovibrational transitions, which implies that the Stark effect due to 
BBR is dynamic. Furthermore, the BBR spectrum covers both the rovibrationally-dominated (long-wavelength) polarizability range and the 
electronically-dominated (short- wavelength) range. This illustrates the need to include both rovibrational and electronic polarizabilities in a 
calculation of dynamic Stark shifts due to BBR. 



from Eq. (5) that this is equivalent to averaging Eq. (7) over 
all M states: 

owH = ^j— Yl a vjM(u), (8) 

leading to a shift AE®f R (T) due to the BBR mean-square 
electric field density (£| BR (u;, T)) of 

i r°° 

A£^(T) = Jf a vJ (u) (^br(w, T)> da;. (9) 

In our model, a v j(uj) involves a summation over terms which 
diverge for frequencies equal to their respective rovibrational 
transition frequencies (Eqs. (2) and (3)). The integration over 
this sum is performed as follows. First, the convergence prop- 
erties of the sum and BBR density function (Eq. (A. 5) in the 
Appendix) allow to interchange the summation and integral 
signs, after which Eq. (9) is evaluated as a series of Cauchy 
principal value integrals. 

We stress that the average polarizability (Eq. (8)) can be ap- 
plied to unpolarized, incoherent electric fields only. To illus- 



trate this, we plot (for v = and J = 1) both the average po- 
larizability a v j(uj) and the polarizabilities for linearly polar- 
ized electric fields «„jm(w) and \M\ = 0, 1 in Fig. 3(a). For 
long wavelengths, cx v jm(u) is dominated by purely rovibra- 
tional transitions. This contrasts the situation for a v j(cj), in 
which the rovibrational contributions to the polarizability av- 
erage out due to the molecular rotation (see also Fig. 4). Sev- 
eral rotational and vibrational transitions occur which decay 
by spontaneous emission (spontaneous lifetime ~ 10 ms 22 ). 
Hyperfine structure (which is ignored in our model) of these 
transitions covers a spectral range of about 1 GHz 13 , which 
would not be visible on the scale of Fig. 3(a). For wavelengths 
shorter than 20 ^m, electronic transitions start to dominate the 
dynamic polarizability, except for narrow spectral regions near 
vibrational transitions where the rovibrational polarizability 
diverges. Another remarkable feature is the absence of cer- 
tain divergences in the J = 1, \M\ = 1 polarizabilities which 
do appear in the J = 1, M = polarizability. This is due to 
the selection rule M' — M = appertaining to electric fields 
linearly polarized along the z-axis (as assumed here). As a 
consequence, states with J = 1, M = are coupled to states 
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Table 4 Dynamic Stark shifts (in mHz) due to T = 300 K BBR for various vibrational states with J = 0, 1 



V 




J = 






J = 1 




Contribution d™j(u)) 


Contribution a^j(u;) 


Total 


Contribution a^j(w) 


Contribution a^j(ai) 


Total 





35 


-27 


8.3 


32 


-27 


4.6 


1 


38 


-33 


5.5 


34 


-32 


1.9 


2 


41 


-39 


1.6 


37 


-39 


-1.9 


3 


43 


-47 


-3.5 


40 


-47 


-7.1 


4 


46 


-57 


-11 


43 


-57 


-14 


5 


49 


-70 


-21 


46 


-70 


-24 


6 


52 


-81 


-29 


49 


-86 


-37 


7 


55 


-111 


-56 


52 


-107 


-55 




(o,o) -<o,i) U aj = i(co) 



200 500 1000 2000 5000 10000 

wavelength / |am 

Fig. 4 (Color online) Long-range wavelength behavior of the 
v — 0, J — 1 polarizabilities shown in Fig. 3(a). Dashed segments 
of each curve correspond to negative-valued polarizabilities, solid 
segments to positive values. Vertical dashed lines indicate the 
position of rovibrational transitions (v, J) — (V, J') coupling to 
J = 0, 2 states. In addition, the polarizability of the (« = 0, J = 0) 
state is shown, which is strictly scalar. Due to the absence of 
rotation, for this state the average polarizability due to rovibrational 
transitions does not average out as for J — 1 states. 

with J' = 0, whereas states with J = 1, \M\ = 1 are not, 
which explains the absence of J' = — J = 1 divergences 
for \M\ = 1 polarizabilities. As expected, the average polar- 
izability a v j(uj) contains all divergences. 

Figure 4 shows the behavior of the dynamic polarizabili- 
ties of J = 0, 1 states at very long wavelengths (electric field 
frequency approaching dc). Here, it is clearly visible that the 
'rotationless' J = state has large polarizability as there is 
no averaging effect by the rotation. In general, the dynamic 
polarizabilities display strong tensorial behavior, in particular 
in cases where the electric field is polarized. This is an im- 
portant feature to bear in mind if Stark shifts due to the radio- 
frequency electric fields used in ion traps are to be considered, 
as these fields have a well-defined polarization. 

3.3.3 Results for BBR shift. As is obvious from 
Fig. 3(b), the Stark effect due to BBR at T = 300 K is dy- 



namic. This situation differs radically from that for atomic 
ions, for which the Stark effect due to BBR radiation can 
be often treated quasi-statically. Thus, the treatment of sys- 
tematic shifts in spectroscopy of HD + must be done with ex- 
tra care, despite the fact that QLS of HD + molecular ions in 
the Lamb-Dicke regime may be done in a similar way as for 
atomic ions 3 ' 7 . 



Dynamic Stark shifts due to T — 300 K BBR to several 
rovibrational levels are calculated by numerical integration of 
Eq. (9) using the Cauchy principal value package of the Math- 
ematica computational program. Results are tabulated in Ta- 
ble 4, in which we also specify the individual rovibrational 
and electronic contributions. The rovibrational contributions 
turn out to produce positive level shifts. This can understood 
qualitatively from Figs. 3(a) and (b). Indeed, the BBR spec- 
trum samples primarily the rovibrationally-dominated spectral 
region (A > 20 /im) where the polarizability attains negative 
values, leading to a positive level shift by virtue of Eq. (2). 
On the other hand, BBR wavelengths below 20 /im primarily 
polarize the electronic structure of the molecule for which the 
polarizability is positive, and which explains the negative shift 
introduced by the electronic contribution (Table 4). We also 
calculate differential BBR shifts to several transitions which 
may be amenable to Lamb-Dicke spectroscopy (Table 5). For 
optical transitions, the differential shifts are relatively small 
and contribute at the level of 10~ 16 . Assuming that the tem- 
perature of the BBR field in an experimental apparatus 5 can 
be determined to within ±10 K, we find from Eq. (9) that the 
BBR shift to optical transitions can be inferred from the polar- 
izabilities derived here with relative accuracy better than 40%, 
or well below 10~ 16 relative to the transition frequency. It 
should be noted that the shifts are much smaller than both the 
HD + hyperfine splittings 23 and Zeeman shifts due to magnetic 
fields typically encountered in experiments 3 . A more refined 
analysis of BBR shifts should therefore include the Zeeman 
effect as well as the hyperfine structure. 
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Table 5 Differential dynamic Stark shifts (mHz) due to BBR at T = 300 K for various rovibrational transitions 



(v',J') 


-(v,J) 


Wavelength/^im 


Contribution aj^(w) 


Contribution a e vJ (oj) 


Total 


Relative /10" 16 


(o,i)- 


(0,0) 


227.98 


-3.9 


0.2 


-3.7 


-28.2 


(1,0)- 


(0,1) 


5.3499 


6.5 


-5.5 


0.9 


0.16 


(4, 1) - 


(0,0) 


1.4040 


7.1 


-30 


-23 


-1.1 


(4,0)- 


(0,1) 


1.4199 


15 


-30 


-15 


-0.72 



4 Conclusion 

The dynamic polarizability of rovibrational states in the 1 s<r s 
electronic state of HD + has been evaluated by combining ex- 
isting data on static polarizabilities with numerical calcula- 
tions done using a simplified model of the HD + molecule. As 
a result of these numerical calculations, new values for radial 
dipole transition matrix elements were obtained which can be 
regarded as an improvement over existing values 8 . The thus 
found dynamic polarizabilities point out that the Stark effect 
due to BBR - an important systematic effect in optical spec- 
troscopy of atomic ions and optical clocks - is highly dynamic 
for the molecular ion HD + , in contrast to BBR shifts to op- 
tical transitions in atomic ions 24 . In this respect, the case 
of HD + is similar to that of neutral molecules 16 . It is fur- 
thermore pointed out that the sign and magnitude of infrared 
dynamic polarizabilities depend strongly on the polarization 
state of the electric fields present. This insight is important for 
the evaluation of another well-known systematic shift in high- 
resolution spectroscopy of trapped ions, namely the Stark shift 
due to the trapping electric fields 25 . Notwithstanding these 
salient features of the HD+ polarizability, it is shown that 
T = 300 K BBR shifts become important for optical spec- 
troscopy of HD+ only at the 10~ 16 level. The smallness of 
the BBR level shifts furthermore suggests that future, more 
refined polarizability calculations should take magnetic-field 
interactions and hyperfine structure into account. 

A Appendix 

Here, we justify the approximations presented in Sec. 3.3.1. 
We start by noting that except the 2pa u electronically excited 
state, all excited-state potential energy curves are located at 
large internuclear range, and that these excited states are con- 
nected to lscr g states by VUV transitions having very poor 
Franck-Condon overlap with lscr s states with low vibrational 
quantum number 17 . Therefore, it is reasonable to assume that 
the larger part of a e vJM {u) stems from bound-free transitions 
from lscrg to 2per u , and that we can use the 2pa u potential en- 
ergy curve of Esry and Sadeghpour 17 to estimate the effect of 
ignoring the dynamic part of the polarizability a% JM {u)). To 
this end, we need to consider the dynamic Stark shift due to 
bound-free transitions. A bound state, subject to an oscillating 



electric field with photon energy E = frcu, will undergo an 
energy shift TiA(E) due to off -resonant bound-free coupling, 
with the corresponding frequency shift being given by 26 

* {m = ±Fvj~^6If. (A.1) 

Here, PV denotes the Cauchy principal value, which is eval- 
uated numerically using the Cauchy principal value package 
of the Mathematica computational program, and T(E)/(2n) 
stands for the bound-free transition rate (in s -1 ) induced by 
an electric field with photon energy E — fko. This transi- 
tion rate can be obtained using Fermi's Golden Rule, an ap- 
proach which was followed by Dunn 27 to calculate cross sec- 
tions a v j(E) for photodissociation of H^. These cross sec- 
tions are proportional to bound-free radial matrix elements of 
the form 



/•OO 




X E f j>(R)D 12 (R) XvJ (R)dR\ 2 , 



(A.2) 

where Xe } j'{R) represents a free (dissociating) state of nu- 
clear motion in 2per„ with asymptotic energy Ef. Ef is re- 
lated to E and the dissociation energy E*j of the bound state 
(v,J)by 

E = E d vJ + E Sl (A3) 

where we have neglected the small (29 cm -1 ) isotopic split- 
ting between the lser 9 and 2pa u dissociation limits 17 . It is 
important to note that the shape of T(E) is governed by these 
wavefunctions via Eq. (A.2), and that the 'dynamic' content 
of the shift A(E) is therefore determined by these wavefunc- 
tions. We calculate XE f j' {R) for the case ofHD + by outward 
numerical integration of Eq. (1) for given energy Ef while us- 
ing the 2p<r u potential of Esry and Sadeghpour 17 . We normal- 
ize the free-particle wavefunctions as done by Dunn 27 , after 
which they may be used to find photodissociation cross sec- 
tions <t v j(E) for various states with v = — 7 and J = 0, 1. 
These cross sections are averages over M levels and therefore 
suited for a treatment of the shift due to BBR (Sec. 3.3.2). 
Multiplying a v j(E) with the flux of photons from the radia- 
tion electric field yields the transition (photodissociation) rate 
T vJ (E) of state (v, J): 

T vJ (E) = 27ra vJ (E)^- = 2na vJ (E)^p-. (A.4) 
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Here, we used the definition of the irradiance I = c€q(£ 2 ). In- 
serting F v j(E) into Eq. (A.l) subsequently produces the level 
shift A vJ (E). 

To test the validity of the approximations made in Sec. 3.3.1 
we apply Eq. (A.l) to two cases. In the first case, we adopt 
the approximation of Sec. 3.3.1 by first deriving the mean- 
square value of the BBR electric field, (£b BR (T)), inserting 
it into Eq. (A. 4), and subsequently calculating the level shift 
in the limit that E -> (i.e. assuming a static field). In the 
second case, we obtain the level shift by proper integration of 
Eq. (A.l) over the BBR energy spectral density. 

For the first case, we find (£bbrCO) fr° m the equation 

^o<£|brCO> = \W{T), 

noting that only half of the integrated BBR energy density, 
W(T), is stored in the electric field. W(T) is found by inte- 
grating the BBR energy spectral density w(u, T)duj: 



W(T) 



- f 

Jo 



w(oj, T)duj 



-du; 



7T 2 (k B T f 

15 (hcf 



(A.5) 



Inserting (^bbr(^)) mt0 Eq- (A. 4), and inserting the resulting 
transition rate T v j(E, T) into Eq. (A.l), we obtain the quasi- 
static approximation to the frequency shift A*£j gg R (T) as 

A^bbrCO= lim .A vJ (E,T). 

Fj — >V 

For the second case, we rewrite the BBR mean-square elec- 
tric field spectral density as 



(£g BR (w,T))dw = -w(uj,T)dto 
eo 

= -J-w(E,T)dE. (A.6) 

After inserting Eq. (A.6) into Eq. (A.4) we obtain the spec- 
trally integrated dynamic BBR shift A^^gp (T) upon evalu- 
ating the expression 

(A.7) 

The errors introduced by the approximation in 
Sec. 3.3.1 can now be simply evaluated from the ratio 
A^ t B C B R (T)/A^ n BBR (T) for various states (v, J) and 



temperatures T. This is possible even so Eq. (A. 2) is incom- 
plete; any numerical prefactor missing there will be common 
to both methods to compute A v j^br{T), and cancel out 
in the ratio. For states with v < 7, we find that the ratio 

I - A^ t BBR(7 1 )/A^ n BBR(7 1 ) < 0.03. Comparing shifts due 
to monochromatic fields in a similar fashion, we observe that 

that the ratio 1 - A vJ (0)/A vJ (E) < 0.1 for A = 4 ^m, and 
decreases to in the static-field limit. This translates directly 
to the accuracy of a v jM{u) claimed in Sec. 3.3.1. 
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